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Abstract. Accretion discs around black holes in X-ray binary stars are warped if the spin axis of the 
black hole is not perpendicular to the binary orbital plane. They can also become eccentric through 
an instability involving a resonance with the binary orbit. Depending on the thickness of the disc 
and the efficiency of dissipative processes, these global deformations may be able to propagate into 
the innermost part of the disc in the form of stationary bending or density waves. We describe the 
solutions in the linear regime and discuss the conditions under which a warp or eccentricity is likely 
to produce significant activity in the inner region, which may include the excitation of quasi-periodic 
oscillations. 
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INTRODUCTION 

Accretion discs involve approximately Keplerian motion around a massive central ob- 
ject, and the general solution allows for nested orbits with smoothly varying inclination 
and eccentricity. The shape of a warped or eccentric disc evolves slowly under the ac- 
tion of stresses in the disc and external forces that deviate from that of a Newtonian point 
mass. Coherent precession of discs is usually possible in binary stars. 

Warped and eccentric discs may be relevant to a wide variety of phenomena in X- 
ray binaries. Long-period modulations such as superhumps and superorbital variability 
may be attributable to the precession of global warping or eccentric modes that are 
not forced but are excited by various instabilities. Sufficiently large discs encounter a 
resonance at the location where the angular velocity of the disc is three times that of 
the binary orbit and may become eccentric as a result [1, 2]. This effect is well known 
in cataclysmic variable stars of mass ratio q < 0.3 [e.g. 3] where, as well as giving 
rise to superhumps, the eccentricity enhances the viscous dissipation and significantly 
affects the outburst dynamics. Similar processes should occur in most X-ray binaries 
with black hole primaries, and indeed superhumps are reported in an increasing number 
of such systems [4, 5, 6, 7, 8]. A warping instability involving radiation forces [9] is more 
likely to occur with neutron star primaries [10] but may also be seen in GRS 1915-1-105, 
which has a very large disc [11]. Stationary warped discs also arise whenever there is a 
misalignment between the orbital angular momentum of the binary and the spin angular 
momentum of the central object [12], a knowledge of which is communicated to the disc 
through general relativistic (gravitomagnetic) or magnetic torques [13, 14]. 

The innermost parts of discs around black holes and neutron stars depart significantly 
from Keplerian motion. Indeed, the rapid relativistic precession of elliptical or inclined 
orbits has often been discussed in connection with quasi-periodic oscillations (QPOs) 



in X-ray binaries, and may best be studied within the context of eccentric or warped 
discs. There is an apparent conflict between the description of the outer part of the disc, 
which supports slowly precessing global deformations, and that of the non-Keplerian 
inner region. In this paper we attempt to make this connection and to describe how a 
global deformation of the disc in the form of a stationary or slowly precessing warp 
or eccentricity may be able to propagate inwards under some conditions, activating the 
inner region and possibly exciting trapped oscillations that may explain high-frequency 
QPOs in accreting black holes. The excitation mechanism itself has been studied by 
Kato [15, 16] and by Ferreira and Ogilvie [17, see also this volume]. 



A small warp or eccentricity can be considered as a perturbation of a standard (circular 
and coplanar) disc. Such a disc supports a variety of wave modes, having a dependence 
on time and azimuth of the form exp(zm0 — i(Ot). In the simplest case of a strictly 
isothermal disc, a local dispersion relation 



can be derived [18], which relates the Doppler-shifted wave frequency <» = (O — mQ. to 

the radial wavenumber k. The integers m and n > are the azimuthal and vertical mode 
numbers, while H is the vertical scaleheight of the disc and Q., K and are the orbital, 
epicyclic and vertical oscillation frequencies characteristic of circular orbits in the given 
potential or metric. All of these quantities depend on the radius R at which the dispersion 
relation is evaluated. The local dispersion relation of a more general disc model can be 
calculated numerically. For this description to be accurate, the wavelength X =2'K/k 
should be much less than R. 

Within this context, a warp corresponds to (m,n) = (1,1) (vertical motion indepen- 
dent of z) and an eccentricity to (m, n) = ( 1 , 0) (horizontal motion independent of z). The 
dispersion relation shows that the warp takes the form of a propagating bending wave 
{k^ > 0) when {o) — Q.)^> max(K"^,^2^) or < min(K"^, while the eccentricity takes 
the form of a propagating density wave when {co — £1)^ > K^. 



A complementary description is provided by theories that consider a warp or eccentricity 
that varies on a length-scale much longer than H and on a time-scale much longer than 
Small-amplitude warps are governed by the equations 
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FIGURE 1. Local wavelength of a stationary warp (left) or eccentricity {right), in units of the local 
vertical scaleheight H, plotted versus the radius in gravitational units, for black holes with spin parameters 
a = 0.1, 0.2, 0.4 and 0.8. The disc is terminated at the marginally stable orbit in each case. 



[e.g. 19], where W describes the amplitude and phase of the inclination of the disc at 
radius R and time t, G refers to a horizontal torque communicated by Reynolds stresses, 
L is the surface density and P = ILH^Q^ is the vertically integrated pressure. These 
equations describe propagating bending waves with essentially the same local dispersion 
relation as in the previous section. (The theories overlap when // -C A -C i?, which 
is possible in a thin disc.) They also allow for viscous (i.e. turbulent) damping of the 
warp, parametrized using a dimensionless number Ow, which is equivalent to the usual 
Shakura-Sunyaev parameter [20] if the disc has an isotropic effective viscosity. 
The simplest equation describing a small eccentricity is 



[e.g. 21], where E describes the amplitude and phase of the eccentricity of the disc 
at radius R and time t, y is the adiabatic exponent and «£■ parametrizes the viscous 
damping. For short wavelengths, this equation also agrees with the local dispersion 
relation of the isothermal disc when 7=1 and Ue = 0. This equation is based on a two- 
dimensional approximation and neglects many of the complications of a shear viscosity 
such as viscous overstability [22]. More sophisticated theories, including nonlinearity 
and all viscous or viscoelastic effects, are available in the literature [23, 24]. 




\J 

dR 



+ R—E + m^{Q}-K^)E (4) 



STATIONARY DEFORMATIONS 



A global warp or eccentricity precesses at only a fraction of the binary orbital frequency, 
so the wave frequency (O is completely negligible compared to Q., k and Q.^ in the inner 
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FIGURE 2. Stationary warp in a disc around a black hole with a = 0.5, e = 0.02 and with aw =0 (top) 
and 0.05 {bottom). Real and imaginary parts of W are plotted as solid and dotted lines. The ampUtude is 
scaled such that W at large r. 



part of the disc. Setting co = 0, we find the local radial wavenumber to be given by 
k^H^ = - K^) {Q? - aj) /a^al in the case of a warp and k^H^ = {Q? - K^) /y^j 
for an eccentricity. Using the expressions for K and D.^ for the Kerr metric [25] and 
taking 7 = 5/3, we find X/H asa function of the dimensionless radius r = Rc^/GM and 
the spin parameter a (see Fig. 1). For a > 0, both W and E propagate at all radii, with 
wavelengths everywhere significantly longer than H. Of all the wave modes described 
by the local dispersion relation, these are the most credible in a turbulent disc because 
of their relatively long wavelengths. 

It is also possible to predict how the amplitudes of the deformations scale with radius. 
In the absence of viscous damping, a WKB analysis of the secular theories shows 
that oc 7?i/^(EH)"i/2 and \E\ oc i?V4(j:^)-i/2 ^pply these results to a steady 
accretion disc in the regime dominated by gas pressure and Thomson opacity, in which 
(assuming a = constant) L oc fV^R-VS and H oc /i/5^2i/20^ ^j^g^.^ / = 1 _ (i?;„/i?)i/2 
[20]. Then \W\ oc f-V^R-i/iO ^nd \E\ oc /-2/5^i/40^ implying a very mild dependence 
of the amplitude on radius. However, the gradients dW /dR and dE/dR do increase 
sharply at small R, because of the rapidly decreasing wavelength. 

When dissipation is taken into account, these solutions are modified by viscous 
attenuation. Our interest is in whether a deformation of the outer part of the disc can 
propagate into the inner region with non-negligible amplitude. Noting that H/R is 
almost independent of R in the above model, we define the constant parameter e by 
H/R = e/^/^r^/^*^. Then the logarithm of the attenuation factor for a complete crossing 
of the disc scales approximately with Uw/s or a^/e. 
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FIGURE 3. Stationary eccentricity in a disc around a black hole with a = 0.5, e = 0.02 and with = 
(top) and 0.05 {bottom). Real and imaginary parts of E are plotted as solid and dotted lines. 



The numerical solutions in Fig. 2 confirm this behaviour. In the absence of dissipation, 
the bending wave reflects perfectly from the stress-free inner boundary and sets up 
a standing wave. Significant attenuation is found when aw exceeds e, in which case 
only a wave with inward group velocity is seen. If aw is increased much further, the 
oscillations are no longer apparent. Only in this case does the inner part of the disc lie 
in the equatorial plane of the black hole as suggested by Bardeen and Petterson [13]. 
Note that W tends to a constant at large r, which corresponds to the inclination of the 
outer part of the disc with respect to the black hole's equator. The oscillatory structure 
of the warp has been noted before [26, 19]. Traces of it may already have been detected 
in numerical simulations [27]. 

The behaviour of the eccentricity is similar (Fig. 3) except for the shorter wavelength 
(still everywhere much longer than H). Again, unless Ue is several times greater than e, 
the eccentricity can reach the inner region. 

Several caveats accompany these solutions. Radiation pressure, which is more impor- 
tant at higher accretion rates, thickens the inner part of the disc, increases the wave- 
length, and reduces the attenuation. Viscous overstability may cause the eccentricity to 
grow, rather than decay, as it propagates inwards. Nonlinearity may be very important. 
The relevant damping coefficients Ow and Ue are not generally equal to a and remain 
relatively poorly understood. In addition, the time for the warp or eccentricity to propa- 
gate into the inner region and to establish the steady solutions shown here can be long. 



CONCLUSIONS 



Accretion discs around black holes in X-ray binary stars may commonly be warped 
or eccentric. Depending on the thickness of the disc and the efficiency of dissipative 
processes, these global deformations may be able to propagate into the innermost part of 
the disc in the form of stationary bending or density waves. This is most likely to occur 
when the disc is hotter and thicker, when the wavelengths are longer and the viscous 
attenuation is less severe. Under these conditions the inner region may be activated and 
trapped oscillations may be excited through nonlinear mode couplings [15, 16, 17]. 
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